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Abstract 

A weak formulation for the so-called semilinear strongly damped wave equation with constraint 
is introduced and a corresponding notion of solution is defined. The main idea in this approach 
consists in the use of duality techniques in Sobolev-Bochner spaces, aimed at providing a suitable 
“relaxation” of the constraint term. A global in time existence result is proved under the natural 
condition that the initial data have finite “physical” energy. 
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1 Introduction 

This paper is devoted to studying the so-called semilinear wave equation with strong damping, namely 

eutt - SAut - Au +f{u) = g, (1.1) 

for e,(5 > 0. The equation is settled in the parabolic cylinder Q = (0,T) x where 11 is a smooth 
bounded domain in V > 1, and T > 0 is a given final time, and is complemented with the initial 
conditions for u and ut and with homogeneous boundary conditions either of Dirichlet or of Neumann 
type. The strong damping is provided by the term —SAup this comes in contrast with the weak 
damping occurring when that term is replaced by +Sut. The function g on the right hand side is 
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a given volume forcing term (here taken of L^-regularity), and the semilinear term f{u) is assumed 
to take the form f(u) = j3{u) — Art, where /3 is a monotone function (more precisely, a monotone 
graph, see Section [2] below) and A > 0. In particular, the internal constraint on u is enforced by 
the non-smooth monotone part /3 of /, whereas the remaining term —\u is related to the (possible) 
nonconvexity of the energy functional associated to the equation. Actually, the main novelty of this 
paper stands in the fact that /3 is assumed to be defined only in a bounded interval Jq of M and to 
diverge at the extrema of /q. A (generalized) function /? with the above properties will be referred to 
as a constraint on the variable u (cf. Section [2] below for more details). It is worth noting that, up to 
purely technical modifications in the proofs, our techniques could be adapted to treat also the case of 
unilateral constraints, i.e., functions /3 whose domain is bounded only from one side. 

Physically speaking, equation (HU appears in a number of different contexts. Let us mention 
here some of them. The main application refers to the study of the motion of viscoelastic materials. 
In this setting, u plays the role of a (scalar) displacement and (11.111 represents the momentum balance 
(where accelerations are included) written in a small strain regime. In particular, respectively in 
space dimensions one and two, the equation describes the transversal vibrations of a homogeneous 
string and the longitudinal vibrations of a homogeneous bar subject to viscous effects. The strong 
damping term —SAut represents the fact that the stress is decomposed in the sum of a pure elastic 
part (proportional to the strain) and a viscous part (proportional to the strain rate), as in a linearized 
Kelvin-Voigt material. We also mention that in the literature, in space dimension three, (HU has 
been introduced to model, e.g., the deviation from the equilibrium configuration of a (homogeneous 
and isotropic) linearly viscoelastic solid with short “rate type” memory (cf. HI] details), in the 
presence of an external displacement-dependent force g — f(u). We do not enter deeper in the modeling 
details, and we refer to ( 22 ] for a physical derivation of models describing the motion of viscoelastic 
media. Let us observe that it would be meaningful to consider here a vectorial (displacement) variable 
u, but we preferred, just for simplicity, to study only the scalar case at least at a first stage. Indeed, 
the extension of our results to the vector-valued case should be possible, at least for constant isotropic 
diffusion, whereas the case of non-constant stiffness (and viscosity) tensors may be somehow more 
involved. In this framework, we also have to quote (possibly adhesive) contact models with unilateral 
constraints (occurring for instance in the case of Signorini conditions) on a part of the boundary. In 
this setting, the (vectorial) operator /3 would force the direction of the trace of u on the boundary in 
such a way to ensure impenetrability (cf., e.g., [H [51 HH])- We are planning to analyze this type of 
models, by using the methods developed in this paper, in future works. 

Equation (HU also appears in the so-called Fremond theory for phase transitions whenever 
microscopic accelerations are taken into account (cf., e.g., [8j|9l[T6]). In that setting, the unknown 
u generally denotes a (scalar) phase parameter, which is related (for a first order phase transition 
in a binary system) to the local proportion of one of the two phases, or components, of a binary 
material. Then, /3 represents an internal constraint forcing u to take values into the physical interval 
whose extrema (often given by — 1 , 1 ) correspond to the pure states, whereas the intermediate values 
represent a mixture of the phases. Physically relevant choices are /3(r) = log(l + r) — log(l — r) 
(i.e., the derivative of the so-called logarithmic potential often appearing in Allen-Cahn or Cahn- 
Hilliard models), or /3(r) = 9/[_i i](r) (i.e., the subdifherential of the indicator function of the interval 
[—1,1], given by /[_ij](r) = 0 for r G [—1,1] and /[_i_i](r) = -boo otherwise). It is also significant 
to consider equation (HU with other kinds of nonlinearities /3 not having the form of a constraint. 
For instance, (HU appears in the recent theory of isothermal viscoelasticity with very rapidly fading 
memory (cf. [12] and references therein), in the sine-Gordon model describing the evolution of the 
current u in a Josephson junction (cf. [25] : there f{u) = sinu), or as a Klein-Gordon-type equation 
occurring in quantum mechanics (then f{u) = \u\'^u for suitable 7 > 0 ). 

Actually, in the case when / is smooth and defined on the whole real line, the mathematical 
literature on equation (HU is very wide (we quote, without any claim of completeness, the papers 
12 HIl m |20l|2ll|24l|2a [32 )■ Referring to [21] for more details, we recall here that one of the 
first essential results on global well-posedness of (the Dirichlet problem for) (11.11) in the 3D case was 
obtained by Webb, who proved in [22 that, if / satisfies standard dissipativity conditions (without 
any growth restriction), then the problem admits a unique strong solution {u, Ut) taking values in the 
space D iJQ(D)] x On the other hand, when one looks for less regular solutions, the 

situation seems different. In particular, it is natural to consider weaker solutions such that the “energy 
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of the system” remains bounded (in the analytical literature this fact corresponds to require that these 
solutions take values in the so-called energy space). Indeed, this type of regularity corresponds to the 
a priori estimate obtained by (formally) testing (11.11) by ut. Then, one can easily realize that, at least 
if the external source g is 0, the functional 

£{u,ut) = -b j(u) - dx, (1.2) 

where j is an antiderivative of f3, tends to decrease in the time evolution. Usually £ is interpreted as 
a physical energy. This is particularly clear in the cases when u represents a displacement (including 
phase-change models where u is related to the effects of displacements at microscopic scales): then 
the component \ut\'^ of the integrand is a density of kinetic energy, whereas the other summands 
correspond to some kind of configurational or potential energy. Consequently, energy solutions can 
be defined as those solutions taking values in the energy space, or, equivalently, keeping finiteness of 
the energy in the course of the evolution. 

From the mathematical point of view, managing this type of solutions may be delicate, espe¬ 
cially in high space dimension, in view of the possibly fast growth of the integrand j{u). Correspond¬ 
ingly, the literature related to this case is much more recent: Kalantarov and Zelik in m consider 
polynomial nonlinearities of the form f{u) ~ without any restriction on the exponent g > 0 

and prove well-posedness of the equation in the energy space. Moreover, they analyze the asymp¬ 
totic behavior of solutions for large values of the time variable. More recently Fata and Zelik in 
have extended these results to the case when f : M. ^ M. is any smooth function satisfying the basic 
dissipativity assumptions, without any restriction on the growth rate (for instance, they may take 
exponential, or even faster growing, terms /). At least up to our knowledge, however, the case when 
/ is of constraint type has never been considered up to now. The typical example we have in mind 
is f{u) = 9/[_i_i](m) — u (cf. (12.41) and (12.51) below), which may describe, for example, some phase 
transition phenomena accounting for microscopic accelerations (cf. [16)1. 

As one addresses the initial-boundary value problem for (HID under our assumptions, the 
main mathematical difficulty comes from the combination of the constraint /3 with the second time 
derivative Uu- Indeed, this feature strongly restricts the available a-priori bounds. To be precise, 
almost all information on the solution has to be extracted from the so-called “energy” estimate, 
i.e., testing the equation by Ut- In addition to that, one can just get some more smoothness of u 
by multiplying (ED by —Au (as is done, e.g., in [IH]). Anyway, this does not help for controlling 
the term uu, which is the main issue from the point of view of regularity. Moreover, the standard 
procedures that one usually adopts for obtaining higher order bounds, like differentiating in time the 
equation, do not seem to work here, at least for a general choice of /3. This seems to be, indeed, the 
main difference of the present problem with respect to first order (in time) equations with constraint, 
for which additional regularity of solutions can be generally deduced by differentiating in time and 
testing the result by ut, whatever is the expression of /3. 

In view of the lack of estimates, we need to build a notion of weak solution which is sufficiently 
general to exist under the sole “energy” regularity. This is, indeed, a somehow delicate issue. In 
particular, one cannot expect to reproduce the same type of results that hold in the case of less 
general nonlinearities /3. To say it shortly, the main novelties of our approach can be summarized in 
two points: 

(i) a relaxed form /?„, of the operator j3 obtained by means of duality techniques; 

(ii) an integrated (both in space and in time) variational formulation where test functions are chosen 
in suitable Sobolev-Bochner spaces. 

These choices permit us, indeed, to prove existence. However, both of them come at some price. 
Namely, 

(i) it will not be possible to intend the equation, and the constraint in particular, in the pointwise 
sense; 

(ii) we cannot exclude the occurrence ot jumps of Ut- Actually, ut may be discontinuous with respect 
to time (and, more precisely, is expected to be discontinuous, as we can show by means of examples). 

However, from a physical point of view, if (ED comes from a variational principle (as the principle 
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of virtual power is), the variational setting in which we introduce the solution is the natural one. In 
particular, the operator (3 (in its relaxed version stands for an internal force which is defined in 
duality with velocities/displaceinents. In addition to that, an internal constraint on the function u is 
still ensured by the definition of the domain of /?„,. Finally, the fact that we can have jumps on the 
velocity Ut w.r.t. time, corresponds to the possible occurrence of internal (or external) shocks, which 
are expected to happen in this framework (cf, e.g., [Ib)b 

A further drawback is concerned with the problem of uniqueness. Actually, we expect the 
occurrence of genuine nonuniqueness, even though some criteria for “physicality” of weak solutions 
may be proposed (cf. Remark 13.31 at the end). 

Let us conclude by giving some more words of explanation for our method. The basic strategy 
of proof is, in a sense, very standard: we replace the singular function /3 by a smooth approximation 
of controlled growth at infinity (e.g., the Yosida approximation), prove existence of a solution u'^ 
to the regularized problem (which basically follows from results already known in the literature, cf., 
e.g., the quoted HZl US [Ml ISH), and then let the approximation parameter e go to 0. Indeed, as a 
consequence of the so-called “energy estimate”, u^, at least for a subsequence, tends to some limit 
u which we would like to identify as a “weak” solution (where, of course, we need to state precisely 
what we mean with this). However, the only uniform bound available for the nonlinear term /3®(it®) is 
in the norm of L^{0,T; and there is evidence coming from concrete examples that we cannot 

go further, at least for general /3. This fact has, indeed, a number of consequences. First of all, 
arguing by comparison, we can obtain an L^(0, T; A')-bound for where Af is a Banach space such 
that L^(fl) and iL^(H) are compactly embedded into X (for example, we can take X = H~^{Vl) 
in the 3D case). This estimate suffices, via a generalized version of the Aubin-Lions compactness 
lemma, to prove strong convergence of in L^(0, T; L^(S7)). However, the limit function ut may 
exhibit jumps with respect to time. Secondly, the limit of can be taken at least in the (weak) 

sense of measures. A crucial point is, as usual, concerned with the identification of its limit. In view 
of our assumptions it looks natural to rely on a suitable version of the so-called Minty’s trick for 
monotone operators, i.e., to combine the weak convergence of u® in some (reflexive) Banach space V, 
the weak convergence of /3®(u®) in the dual space V', and a lim sup-inequality. A look at the estimates 
suggests that an admissible choice for this procedure is (in the Dirichlet case) the Sobolev-Bochner 
space V = iL^(0,T;L^(H)) n L^(0, T; FIq (H)) (in the Neumann case, Hq{Q) is simply replaced by 
iJ^(H) and no further difficulties arise). In such a setting, the constraint (3 has to be reinterpreted in 
a relaxed form Pw acting as a maximal monotone operator from V to 2^ (cf. Definition 12.131 below; 
see, e.g., [TTJ|TS] for some additional background). Correspondingly, equation (II. 1|) has to be intended 
as a relation in V'. Let us point out that, from a physical point of view, in the case when 
corresponds to a mechanical balance equation (i.e., to the momentum balance equation), our weak 
formulation takes the meaning of a duality between forces and velocities in time and space (see [7] for 
a similar approach, but in a different setting). To avoid occurrence of second time derivatives in the 
weak formulation, we also need to integrate by parts with respect to time the second order term utt 
(cf. (I2.31|1 L Actually, these modifications will permit us to solve our original problem on the whole 
time interval (0,T), but also to write “pointwise” the duality relation in any subinterval (0,t), with 
a physically consistent interpretation of the corresponding constraint. Finally, an energy inequality is 
proved to hold on (almost) every subinterval of [0,r]. We end observing that the behavior of weak 
solutions (at least in the homogeneous Neumann case) may be clarified by considering a spatially 
homogeneous setting. For instance, if f{u) = and g = 0, (11.11) reduces to the prototype 

ODE Utt + i9/[_iq](u) 9 0 whose solutions can be easily described, especially in relation with the 
jumps of Ut (cf. Remark 1 2. 41 for more details). 

The remainder of the paper is organized as follows: in Section|2]we introduce some amount of 
preliminary material mainly related to maximal monotone operators and duality methods; moreover 
we present the notion of weak solution and state the related existence result. Then, the proof is 
detailed in Section [31 where we also give a number of remarks illustrating our results at the light of 
simple finite-dimensional examples. 
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2 Preliminary notions and main result 

Let C be a smooth bounded domain (with TV > 1) of boundary L and let us consider the 
interval [0,T], for some fixed hnal time T > 0. Let us set H := Lp'iyi) and use the notation (•, •) for 
the scalar product both in H and in . Let also the symbol || • || denote the corresponding norms. 
In our analysis, we will consider either Dirichlet or Neumann boundary conditions for (HI); hence we 
introduce a notation suitable for addressing both cases in a unified way. So, we put V := H^{n) in 
the Neumann case, and V := in the Dirichlet case. In both cases, V will be endowed with the 

standard (Sobolev) norm, indicated by || • ||v. Moreover, we will denote by (•, •) the duality pairing 
between Y' and V. In general, we will indicate by || ■ ||x the norm in some Banach space X (or in 
X^). 

We let A stand for the weak form of (minus) the Laplace operator seen as an unbounded linear 
operator on H whose domain D{A) depends on the boundary conditions. Namely, in the Neumann 
case, we set 

Av:=-Xv, D{A) := Hlin), (2.1) 

where denotes the space of the iL^-functions having zero normal derivative (in the sense of 

traces) on d^l. Correspondingly, in the Dirichlet case, we set 

Av:=-Av, D{A) := H‘^{n)nV = (2.2) 

In both cases, ^ is a positive operator (strictly positive for Dirichlet conditions). Morever, A can be 
extended to the space V by setting (for both choices of boundary conditions) 

{Av,z)= / Xv - Xz dx. (2.3) 

Jn 

This extension, which turns out to be linear and bounded from V to V', will be identically noted as 
A; indeed, we believe that no danger of confusion exists at this stage. 

Next, we specify our assumptions on the semilinear term f{u). First, we suppose that / may 
be decomposed as 

f{u) = I3{u) - Xu, (2.4) 

where A > 0 and /3 is a maximal monotone graph in R x K such that 

0€/3(0). (2.5) 

Indeed, just for simplicity and with no loss of generality, we require the closure of the domain of /3 to 
be the interval [—1,1]. In addition, it is not restrictive to assume the normalization 0 € /?(0), which 
turns out to be useful especially in the Dirichlet case. 

Referring the reader to mm for a complete survey on the theory of maximal monotone 
operators in Banach and Hilbert spaces, we just observe here that, thanks to (12.5L there exists a 
convex and lower semicontinuous function j : —> [0,+oo] such that /? = dj, D{j) = [—1,1], and 

j(0) = minj = 0. Here, D{j) denotes the domain of the convex function j, i.e., the set where j takes 
finite values. 

It is well known that the graph (3 induces maximal monotone operators (identically noted as 
/3 for simplicity) both in H and in L^{Q), where Q := (0,T) x D. For instance, one has ^ S P{u) in 
the H-sense if and only li u,^ G H and ^{x) £ I3{u{x)) for a.e. a; G D. Moreover, let us define the 
convex functional 

J : H ^ [0,+oo], J{u) := f j{u) dx, (2-6) 

Jo. 

where the integral may well be +oo in the case when j{u) ^ L^(D) (i.e., when u ^ D{J)). Then, 
P = dJ in H, namely the operator induced by /I on iL coincides with the iJ-subdifferential of the 
convex functional J. As is customary when dealing with multivalued operators, we shall often iden¬ 
tify maximal monotone operators with their graphs (cf., e.g., [31 HH])- With the above notation, 
equation (EH), where the coefficients e and 5 have been set to 1 for simplicity, becomes 

Utt + Aut + Au + P{u) — Xu 5 g in (0, T) x Q. (2.7) 
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Note the occurrence of the inclusion sign, motivated by the fact that (3 may be multi-valued. 

In view of (EH) (or of (EH), (EE) can be read as a relation holding in L^(0, T; i7) (and thus 
interpreted as a pointwise inclusion almost everywhere in Q). Indeed, (12.71) looks as the most natural 
and appropriate weak formulation of the strongly damped wave equation in the case when /3 is a 
smooth monotone function defined on the whole real line. On the other hand, though ( 1 ^ is still 
perfectly meaningful from the mathematical viewpoint under our assumptions E1-E3, proving 
existence of solution in the current setting seems to be out of reach (see Remark 12.41 below for a 
counterexample in the spatially homogeneous case). Mainly, what seems to fail is the possibility to 
interpret point-by-point the equation, and in particular the constraint j3. 

Hence, we need to construct a furtherly relaxed formulation of the equation, for which one 
might be able to get existence. In performing this program, we would like our new concept of solution to 
be still somehow physically consistent. Namely, weak solutions should comply with thermodynamical 
principles (like the energy inequality), satisfy a proper form of the constraint, and be obtained as limit 
points of families of functions solving physically sound regularizations of the equation. To start with 
this program, we set 

V := i7^(0,T;iJ) nL2(0,T;H). (2.8) 

Note that, in view of standard results on vector-valued functions, the above space coincides with 
H^{Q) in the Neumann case. The duality pairing between V' and V will be noted by ((•, •)). We also 
consider the space "H := L/^iQ) = T^(0, T; H) endowed with the natural scalar product, noted here as 
((•, •)). Thanks to standard results on Sobolev spaces, the inclusions V C 'H C V hold continuously 
and densely provided T-i is identified with its dual by means of the above scalar product. Actually, 
the weak formulation of our problem will strongly rely on the parabolic Hilbert triplet V). 

We also need similar concepts in the case when the time interval (0, T) is replaced by (0, t) for 
0 < t < T. Namely, we set Qt ■= (0, t) x H, and, correspondingly, we note by {{v, 2 ))(o,t) the (standard) 
scalar product in Jit ■= L^{Qt) and by ((•,-))(o,t) the duality between Vt := (~lL^( 0 ,t;H) 

and its dual. We also set 

Vt,o := {u e Vt : u = 0 on {t} x O}, (2.9) 

where relation u = 0 is intended in the sense of traces (in time). Clearly, Vt,o is a closed subspace of 
Vf. Then, it ip G Vt,Q, extending it by 0 for times larger than t, we obtain an element of V, noted in 
the following as ip. Correspondingly, it p G V', we can naturally define its restriction p* to the time 
interval ( 0 ,t) by setting, for ip G Vt,o, 

{{v*,‘P))io,t) ■■= ( 2 . 10 ) 

Actually, it is readily checked that p^ G V( q. Moreover, the restriction operator p i-G- p* is linear and 
continuous from V' to q. 

With the above notation at disposal, we extend the functional J to time-dependent functions 
by setting (see (12.61) 1 

J : 7i ^ [0,+oo], J{u) := f J{u)dt= f f j{u) dx dt, (2-11) 

Jo Jo Jfl 

where, as before, the integral may also take the value -l-oo. Analogously, for t G (0,T], we put 

J{u) ds = f f j{u) dx ds. (2.12) 

Jo Jn 

As noted above, the "H-subdifferential dj{u) (or the analogue for J{t){u)) can be still interpreted in 
the “pointwise” form f3{u). 

We are now ready to introduce the weak form of the constraint /3. We shall present most 
of the construction by working on the time interval (0,r). The adaptation to subintervals (0,f) 
is straighforward and we mostly leave it to the reader because we do not want to overburden the 
notation. That said, we start by setting Jy := J\v- It is readily proved that Jy is convex and 
lower semicontinuous on V. Hence, we may take its subdifferential with respect to the duality pairing 
between V and V'. Namely, for ^ gV and u £ V, we put 


•Jp) '■ J~Lt [0,-boo], f 

Jo 


i G Pw{u) u GV, C G 12', and {{^,v — u)) + Jv{u) < Jy(u) for all v GV. (2-13) 
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The idea of “relaxing” /? in this way is not new; for instance, the same method has been applied 
in la na [30] in other contexts. It is worth noting from the very beginning that u S D{/3w) still 
implies u £ [~ 1 , 1 ] almost everywhere; in other words, the weak operator still forces u to assume 
only “physically meaningful” values. Note that an alternative, but essentially equivalent, approach 
based on variational inequalities has been devised in m for the Cahn-Hilliard equation with dynamic 
boundary conditions. The novelty occurring in our case is related to the use of “parabolic” (Sobolev- 
Bochner) spaces. Indeed, this choice seems particularly appropriate for the present problem as far as 
it permits us to overcome some issues related with the (expected) low regularity of weak solutions. 

Let us now characterize a bit more precisely the operator We follow here the lines of 
[TT1[TH] (see also [5]). Firstly, we observe that (see, e.g., [B] Prop. 2.3]), if u G V, ^ £ H, and ^ £ /3(it) 
a.e. in Q, then ^ £ (3w{u)- Namely, if /3|v denotes the restriction to V of the “pointwise” operator /3, 
then Pui extends /3|v. In other words, the “strong” constraint implies the “weak” one. Moreover 
(cf. [3 Prop. 2.5]), 


if M £ V and ^ £ (3w{u) fl then ^ £ (3{u) a.e. in Q. 


(2.14) 


In general, however, the elements ^ £ Pw{u) (which lie, by definition, in the space V') need not belong 
to H. Hence, the graph inclusion /3jv C Pw is generally a proper one. Nevertheless, if ^ £ then 

^ “automatically” gains some more regularity. 

In order to explain this phenomenon, we proceed along the lines of [SO] Sec. 2]. Namely, for 
t £ (0, T], we set 

Xt := C^iQt), for Neumann boundary conditions, (2.15a) 

Xt := {u £ C^{Qt) : m = 0 on [0,t] x P}, for Dirichlet boundary conditions. (2.15b) 

For t = T we simply write X = Xt- We also set, in both cases, 

<^t,o ■= {v G Xt : = 0 on {t} x O}. (2-16) 

The space Xt (hence its closed subspace Xt^o) is naturally endowed with the supremum norm jj • jjoo. 
Moreover, also thanks to the smoothness of ft in the Neumann case, Xt fl Vt is dense both in Xt and 
in Vt- Let now ^ £ V' (the analogue applies with straighforward modifications to ^ £ V() and let us 
suppose that if restricted to the functions (^ £ n V, is continuous with respect to the <Y-norm, 
i.e., there exists C > 0 such that 

-z))| < C'lkiloo for any 2 ; £ df n V. (2-17) 


In that case, by density, ^ extends in a unique way to a bounded linear functional on X. Namely, there 
exists a unique T G X', which can be seen as a Borel measure on Q in view of Riesz’ representation 
theorem, such that 




dr 


for any z G X nV. 


(2.18) 


In this situation we say that the measure 7~ represents ^ on X. Actually this situation automatically 
occurs when ^ is an element of a weak constraint. Indeed, by an easy adaptation of [301 Prop. 2.1] 
(which, in turn, is based on the results of [H]), one can see that, up to some adjustment related to 
the boundary conditions, any ^ £ Pw{u), when restricted to continuous functions, is represented by 
a measure T defined on the parabolic cylinder Q. Such a measure, in turn, is related to the original 
operator j3 in the following way (cf. [TTl Thm. 3] for further details): noting as T = Ta + Ts the 
Radon-Nikodym decomposition of T, with Ta {Ts , respectively) standing for the absolutely continuous 
(singular, respectively) part, we then have 


TaUGL\Q), (2.19) 

Ta{t,x) G P{u{t,x)) for a.e. {t,x) G Q, (2.20) 

{{?:«))- [ [ TaU dxdt = sup \ [ [_z dTs, 2 : £ A, z((3) C [-1, 1]1. (2.21) 

Jo Jfi {JJq J 
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Hence, the continuous part Ta of the measure T satisfies the constraint pointwise (in view of (12.201) '). 
whereas the singular part Ts is characterized by (12.211) . 

In particular, we expect that condition (12.211) could be made more precise. Namely, noting as 
Ts = p\Ts\ the polar decomposition of 7^, where \Ts\ is the total variation of Ts-, proceeding along the 
lines of m Thm. 3] one may prove that 


p e |7^|-a.e. in Q. 


( 2 . 22 ) 


In other words, we expect the singular part of T to be supported on the set where |u| = 1 and that 
p = 1 where u = 1, p = —1 where u = —I. In this sense, also the singular part of T is, at least 
partially, reminiscent of the expression of the graph /3. 

Actually, the characterization (12.221) is proved in [T5] in the case when V = i7o(H), H a 
bounded domain of and may be likely extended to the present situation. However, a detailed 
proof may involve some technicalities particularly related to the facts that we are working in the 
parabolic cylinder and should distinguish between the Dirichlet and Neumann cases. For this reason, 
we omit details here. We note, however, that (I2.22p is straighforward whenever we additionally know 
that u G A n V (i.e., u, beyond lying in V, is continuous). Indeed, in that case, from (12.211) there 
follows 



pu d|7;| 



U,u)) 


TaU dx dt 


lo JQ. 


= sup 



z d7^, z & X, z{Q) C [-1,1] 


ir.i(g), 


(2.23) 


the latter term denoting the total variation of the measure T- Comparing terms, we then deduce 
pu = 1 |7^|-a.e. in Q, as desired. 

Remark 2.1. It is worth observing that, in the Neumann case, the singular component T of the 
measure T representing ^ G Pw{u) may be, at least partially, supported on the boundary of Q. Let us 
see this by a simple one-dimensional example. Let Q = (—1,1), uix) = x, j = Then, if ^ = ol5i 

for some a > 0, where is the Dirac delta concentrated in 1, it is clear that, for any v G 1) 

such that V G D{J) (i.e., such that —1 < v{x) < 1 for all x G Q), there holds 


(^, V — u) = Qf(u(l) — u(l)) = a(v(l) — 1) < 0 = J{v) — J{u). (2.24) 


Hence, ^ G l3-w{u) by definition of subdifferential. 

Having clarified the nature of the weak constraint /3^, we can now observe that equation (EID admits 
a natural energy functional 

£{u,ut) = -b j(u) - da: = i(||utf + \\S/uf + 2J{u) - A||u|p). (2.25) 

Indeed, testing (EH) by Ut and integrating in time and space, one can get that the value of £ at any 
time t > 0 is bounded by the initial value £{uq, ui) plus the power of the external applied forces g (see 
(|2.41l) below). Actually, as will be explained later on, the low regularity of solutions does not allow 
us to perform this estimate directly for weak solutions, but only for a suitable approximation of the 
problem. This is basically the reason for which we will only be able to prove an energy inequality for 
weak solutions, cf. Theorem 12.51 below. 

We can now introduce our assumption on the source term and on the initial data, the latter 
corresponding exactly to the finiteness of the “initial energy”: 

gGL^iO,T-H), (2.26) 

uq G H, ui G H, J{uo) < oo. (2.27) 

Then, we can make precise our concept of weak solution (to be precise, we shall speak of “parabolic 
duality weak solution” or something similar, but we will rather use “weak solution“ just for simplicity): 
















Definition 2.2. A couple (m,^) is called a weak solution to the initial-boundary value problem for 
the strongly damped wave equation with constraint whenever the following conditions hold: 

(a) There hold the regularity properties 

ueW^’°^{0,T-H)nH\0,T;V), (2.28) 

^GV'nA"; (2.29) 

more precisely, by (12.2911 we may require that there exists a measure T € X' representing ^ over 
X in the sense of (12.181) . 

(bl) There holds the following weak version of (12.71) : 

- {{ut, ifit)) + {ut{T),^{T)) + + {{Vu, Vip)) 

+ - Kiu,p)) = {ui,ipiO))(ig,ip)), V(^eV. (2.30) 

(b2) An analogue of ()2.30l) holds also on subintervals, in the following sense: for any t G (0, T] there 
exists a functional ^(t) G V[ such that 

- {{ut,'Pt)){o,t) + {ut{t),'p{t)) + ((Vut, V:^))(o,t) + ((Vu, V(/:))(o,t) 

+ ((C(i)’ 'i^))(o,t) “ T))(o,t) = (■“!: ¥^(0)) + ((g, <p))(o,t)j Vcp e Vt- (2.31) 

Moreover, lies in V[ D if/; hence it is (uniquely) represented over Xt by a measure Tp). In 
addition to that, for every t G (0,T], the functionals ^ and are compatible, namely, for every 
ip G have (p denoting the trivial extension of p) 

Uw,T))(o,t) = U,T))- (2.32) 

In other words, the functional when computed on the elements ofVty, coincides with the 
canonical restriction of ^ (cf. (12.1011 ). 

(c) There holds the inclusion 

C G /3^(m) in V. (2.33) 

More generally for every t G (0,T], G llw,{t){u) in V/. Here represents the weak version 
of P in the interval (0, t); namely, is the subdifferential of the restriction of to Vt with 

respect to the duality product between Vt and V/. 

(d) There holds the Cauchy condition 

u\t=o = uq. (2.34) 

(e) For every s,t G [0,r], the couple (u,^) satisfies the equality 

~ \Wt\\'L:‘{s,t-H) + {nt(t),u{t)) + + \\^'>^\\‘L2(^s,t-H) + ((?(t))''^))(0,t) “ {{C(s), u)) (^o,s) 

- = (Ms), u{s)) + ^I|Vm(s)|P + ^ X 

It is worth discussing a bit how the above formulation has been obtained from (EH). First of all, /3 has 
been replaced with its “relaxed” form . Correspondingly, (E3) has been restated in the “parabolic” 
dual space V' by using the test function p G V and performing suitable integrations by parts. In 
particular, a key point stands in the integration in time of the “hyperbolic” term utt- Indeed, no 
second time derivatives of u appear in (12.301) (or in (12.311) '). In addition to that, the Cauchy condition 
for Ut is now “embedded” into (12.301) and (12.311) . 

Remark 2.3. We need to explain in some detail the “meaning” of (12.3111 . especially in relation with 
the constraint term. Actually, if ^ G there is no canonical way of restricting ^ to obtain an element 
of V/. The best we can do is restricting ^ as explained in (12.101) to obtain a functional G V/ q. 
However, writing ()2.31l) as a relation in V/q (i.e., considering only test functions p G Vt,o) would 
give rise to some information loss. Namely, it may happen that the singular part Ts of T is, at least 
partially, supported on some set of the form {t} x H (or, correspondingly, {{^(t), T))io,t) may also depend 
on the trace of p on {t} x H). 


9 
























Remark 2.4. It is worth noting that, according to the above definition, ut need not be continuous 
with respect to time, independently of the target topology. This fact is a distinctive feature of this 
problem and there seems to be no hope of avoiding jumps of Ut, at least for a general constraint /3. 
Here is a simple example where a jump occurs. Let us consider the case of spatially homogeneous 
solutions to the Neumann problem with A = 0 and g = 0. In other words, we reduce our problem to 
the “toy model” represented by the ODE 


w«+/3(m)9 0, (2.36) 

a weak solution u to which exists according to our theory. Let us also choose /3 = Then, if 

we take, for instance, uo = 0 and ui = I, we get that u{t) = t at least for t £ [0,1). As t gets to 1, Ut 
must develop a discontinuity, otherwise, u(t) would become strictly larger than I for t > 1, and the 
equation would no longer make sense. Hence, the only possibility for the trajectory u{t) is to jump 
instantaneously in such a way that, in a right neighbourhood (1,1 + e) of t = 1, 

ut(t)=i, u{t) = l + i{t-l), £<0. (2.37) 


The trajectory, at least in principle, may “choose” at which level £ the time derivative Ut “decides” to 
jump (hence we have no uniqueness). If it jumps to £ < 0, then u{t) starts to decrease from the value 
I at a constant velocity £ until it reaches the value —1 (where a new jump of Ut must occur). On the 
contrary, if ut jumps to £ = 0, then it will be either u(t) = 1 and Ut{t) = 0 forever, or after some 
time Ut may make a further jump to some £ < 0, starting from which u begins to decrease as specified 
above. More precisely, we can notice that, for (12.361) . the weak formulation over (0,T) (cf. ()2.301) 1 
reads 

T 

umdt + g>{T)ut{T) + {{^,g,))=g,{Q)ui V^£H\0,T). (2.38) 



Hence, it is easy to check that, for all £ < 0, the function u described above solves (I2.38P on a suitable 
interval [0,r] with T > I chosen sufficiently small so that no other jumps of Ut occur. Note in 
particular that different choices of £ correspond to different “values” of ^ G fiw{u). Indeed, from (12.381) 


we get 

pT 



- (fitdt- £v?t d£ + £<p(T) + ((^, (/?)) = <p(0). 

(2.39) 


Jo Ji 


whence 

U,:p)) = {l-£M1), 

(2.40) 


or, in other words, ^ = (1 — £)5t=i (<5 standing for the Dirac delta) and we can notice that this is 
consistent with the above characterization of Indeed, at least for £ < 0, £ = 1 is the only time at 
which u takes the value 1 and ^ may have, and in fact has, a “singular” part. However, we will see 
in the sequel (cf. Remarks 13.71 and 13.81 below 1 that not every jump of ut (or, in the current example, 
every value of £ > 0) is “physically” admissible. 


We can now introduce the statement of our main result: 


Theorem 2.5. Let us assume (12.41) . (12.51) . (12.261) . and (12.271) . Then, there exists at least one weak 
solution (m, ■f), in the sense of Definition 12.21 to the initial-boundary value problem for the strongly 
damped wave equation with constraint. Moreover ut £ BV{0,T; X) for any Banach space X such 
that L^{Ll) and V' are compactly embedded in X. 

In addition, for almost every s £ [0,T) (surely including s = 0) and every t £ (s,T], the following 
version of the energy inequality holds: 

£(u(t),ut{t)) + < £:(M(s),Mt(s)) + f {g,ut) dr, (2.41) 

J S 

where £ is defined in (12.251) . 

Finally, in the case when we additionally have 


uo £ D{A), 


(2.42) 
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then u enjoys the additional regularity property 


ueC^{[0,T]-D{A)). (2.43) 

Namely, u{t) belongs to D{A) for every t £ [0, T] and 1u{t) is continuous when the target space is 
endowed with the weak topology. 

3 Proof of Theorem 12.5 

step 1. Approximation. We start by introducing a natural regularization of (the strong form of) 
equation (Ell) depending on an approximation parameter e (which will then be let go to 0). To this 
aim, for e £ (0,1), we let : R ^ [0,oo) denote the Moreau-Yosida regularization of j (cf., e.g., 
nni for details). In particular, j® turns out to be convex and lower semicontinuous. Moreover, its 
derivative [3'^ := dj^ corresponds to the Yosida approximation of /3 = dj. Under our assumptions /3^ 
is monotone and globally Lipschitz continuous on the whole real line and it satisfies /3^(0) = 0. We 
also set ^ 

J^{u) := / j^{u) dx, J^{u) ■.= f f dx dt. (3.1) 

Jfi Jo Jn 

Moreover, we regularize the initial data by taking, for e £ (0,1), ug and u\ satisfying 

ul&D{A), u\&V, J^{uq) < J{uo), (3.2) 

Uq —>■ uq in V, ul —)• ui in H. (3-3) 

The construction of approximate initial data complying with (I3.2I) - (I3.3D is standard. For instance, one 
may take Uq as the solution to the elliptic singular perturbation problem 

Uq £ D{A), Uq + eAuq = Uq in H. (3.4) 

In particular, the last of (13.21) can be shown by testing the equation in (13.41) by /3'^(ug) and noting 
that 

(P^Uq), Uq - Uq) > J%Uq) - J^{uq) > J%Uq) - J{uq), (3.5) 

the latter inequality following from the monotonicity of the Moreau-Yosida regularization with 
respect to e. 

We are now ready to introduce our approximated equation 

ul^ + Aul + Au’^ + {u^) — \u’^ = g in (0,T) x U. (3.6) 

Correspondingly, we have the following well-posedness and regularity result: 

Theorem 3.1. Let us assume (12.41) . (12.51) . (I2.26|) . and (12.271) . For e £ (0,1), let Uq, uf, and (3^ he as 
detailed above. Then, there exists a unique solution 

u^ £ If^(0,T;If)nW^’°°(0,T;V)nLr^(0,T;D(A)) (3.7) 

to equation dMl), complemented with the initial conditions u^\t=o = Uq and u^\t=o = uf. Moreover, 
for all t,s G [0, T], the following energy equality holds: 

bumr + [ f{u%t)) dx - ^\\u%t)r + hNu^t)r + 

J Q 

= ^ll«K5)f + [ /(«^(,))dx-^||u^(s)|p + i||Vu-(s)f + [\g,ut)dr. (3.8) 

The proof of Theorem 13.11 is fairly standard and could be carried out, e.g., by following the lines 
of [23]. Here it is just worth noting that the regularity conditions stated in (13.71) are compatible 
with the assumptions (I3.2D - (I3.3D on the regularized initial data. Moreover, one could easily check 
that (E3 can be (at least formally) obtained testing (13.6|) by -|- Au^ and performing integrations 
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by parts. In particular, the term can be managed thanks to the Lipschitz continuity of Yosida 

approximations. It is also worth noting that, in this regularity setting, equation (EH) makes sense 
pointwise; indeed, all its single terms belong to the space (in particular, we do not need 

to regularize the source term g: condition (12.261) is enough). Hence, testing the equation by uf is 
allowed: this gives relation (13. 8p by means of well-known chain rule formulas. 

Step 2. A priori estimates. We now derive a number of bounds, uniform with respect to the 
regularization parameter e, for the solutions given by Theorem 13.11 First of all, testing (1^ by 
S Vt, integrating over Qt, and performing suitable integrations by parts (both in space and in time), 
we deduce the integrated (weak) formulation 

- {{uh‘Pt))io,t) + -f ((Vu)^, V<p))(o,t) + ((VM^ Vv9))(o,t) 

+ ((/3^(?^''),</?))(o.t) - A((u^v?))(o.t) = (ui,v?(0)) + ((g,</?))(o.t), e Vt. (3.9) 

Of course, (13.91) holds in particular for t = T and v? € V. Next, setting s = 0 in (13.81) . or, in other 
words, testing EH by uf and integrating over Qt, t S (0,T], we find (by Young’s inequality) 

\\\nmf + [ dx - + i||Vu^(t)f + = M,{e) + 

J O 

< Mi{e) + \\g\\\' 2 (^Q^f,H) + lli2(o,t;ff) =■ |li 2 (o,t;i/)i (3.10) 

where we have set 

Mi[e) := \\\u\f + J-(ug) - + i||Vugf, Nh[e) := M,[e) + \\g\\l.(,,,H). (3.11) 

and we may notice that actually Mi and M 2 are bounded uniformly in e due to (13.2^ - 03.3L (12.271) . 
and (I2.26L Let us also observe that, thanks to the properties of the Yosida approximation (cf. [10)1. 
there exists a constant c > 0 such that j®(r) — Xr^ > —c for all r € K. and all e G (0,1). Hence, 
applying Gronwall’s lemma to (I3.10p . we obtain 

|K(t)|| < M for all t G [0,r], (3.12a) 

and for all e G (0,1). Here and below M denotes a positive constant, possibily different from line to 
line, depending on the problem data, but independent of e. From (13.1011 we also get 

lk^llffi(o.T;y) < (3.12b) 

[ dx < M for all t G [0,T], (3.12c) 

Jn 

for all e G (0,1). Setting now (p = u'^ in (13.911 and taking t = T, we deduce 

/ / I3^{u^)u^ dx dt-k ||Vu^||i 2 ( 0 .T;ff) + ^||VM''(T)f = ||uf||i 2 ( 0 ,T;H) “ {ul{T),u%T)) 

J0 Jq ^ 

+ + + {{g,u^)). (3.13) 

Now, the right hand side is bounded uniformly in e due to (13.1211 . (I2.26I1 - (I2.27I1 and (I3.2|1 - (I3.3L 
Moreover, it is easy to check (cf. also [271 Appendix]) that there exist constants ci > 0, C 2 > 0 
independent of e G (0,1) such that Ci|/3®(r)| < f3^{r)r -|- C 2 for all r G M. Hence, (13.1311 entails 

ll/3^(w^)l|Li(o.T;Li(n)) < M, (3.14a) 

for all e G (0,1). Then, using once more (13.1211 and comparing terms in (13.611 . we also find 

\\uu\\l^{o,t-x) < M, (3.14b) 
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for any Banach space X such that L^(r2) C X and and V' C X with continuous and compact 
embeddings. In particular, it is not restrictive to assume X be the dual of a reflexive and separable 
space (for instance, in dimension iV = 3, one may take X = Hence, we have 

Ikt llvyi.i(o.T;X) < M, (3.14c) 

for all e € ( 0 , 1 ). 

Step 3. Passage to the limit. Now we aim at letting e \ 0. From (l3.12aD - p.l2b|) we deduce that 
there exists a function u of the regularity specified in (12.281) such that 

^ u weakly in iJ^(0, T; H), (3.15a) 

^ u weakly star in {0,T; H), (3.15b) 

and in particular 

u^{t) u{t) strongly in H for all t G [0,T], (3.15c) 

u^{t) u{t) weakly in V for all t G [0,T]. (3.15d) 

It is worth stressing that the above convergence relations, as well as the ones that will follow, are 
intended to hold up to extraction of suitable (nonrelabelled) subsequences of e \ 0 . 

Since V is compactly embedded into i/, in view of conditions (I3.12bl) and (|3.14cl) we can apply 
[HI Corollary 4] with the three spaces V CC H C X and p = 2 in order to obtain that 

Uf ^ Ut strongly in L'^{0,T\H). (3.15e) 

Moreover, condition (l3.14cD implies that the functions rtf are uniformly bounded in BV{0,T; X) (for 
the properties of vector-valued HH-spaces one can refer, e.g., to [TUI Appendix]). In view of the fact 
that we may assume X be the dual of a reflexive and separable space, we can employ a generalization 

of Helly’s theorem (cf., e.g., Thm. 3.1] or [TUI Lemma 7.2]), providing a function v G BV (0, T ; X) 

such that 


u\{t) v{t) weakly star in X for all t G [0,r]. (3.15f) 

It is easily seen that v coincides with ut almost everywhere. Hence, up to changing the representative 
of Ut, we may assume v = Ut everywhere on [0,r]. Moreover, combining (I3.15fp with (|3.12aL we 
obtain 


Moreover we get 


u^{t) Ut{t) weakly in H for all t G [0,T]. 


Ut ^ Ut weakly star in BV{0,T;X). 


(3.15g) 

(3.15h) 


Let us now show that the functions f5^{u^) are uniformly bounded (with respect to e) in V'. Actually, 
writing p.9l) for t = T, and using Holder’s inequality, (13.3p . (|2.26l) and the estimates (I3.12al) and 
p.l2b|) . we find 

|((^^(u^), (p))| < ||Ut ||L2(o,T;H)||‘/5t||L2(o_T;ff) + Ikf (F) ||||(p(T) || 

+ II 11^2(0,T;H) II V(p||i2(o,T;ff) + II Vm® ||^2) || Vp||^2) 

+ ^l|w^lU2(o,T;ff)||<PlU2(o,T;H) + ll^f || ||<p(0) || -|- Hff || L2 (o,T;H) || V^ll (0,T;ff) < C'II'PIIv, (3.15i) 

for all ip G V, where C > 0 is independent of e. Therefore, we can infer that there exists € V' such 
that 

/3'^(ti^) ^ ^ weakly in V'. (3-15j) 

Next, recalling the definition (12.1511 of A, from (I3.14al) we obtain that there exists a measure T G X' 
such that 

f3'^{u'^) T weakly star in X'. (3.15k) 
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In view of the density oi X HV both in X and in V, the measure T represents ^ on X, i.e. (12.1811 holds. 

Let us now go back to dSH), now rewritten for general t G (0, T] and (p G Vf Then, rearranging 
terms, and using the above convergence relations (I3.15p . we obtain that there exists the limit 




lim 
e\,0 L 


(K,<Pt))(0.t) - K(i),‘p(i)) - ((Vu^,Vip))(o.t) - ((VM®,V<p))(o,t) 


+ M(^^7‘p))(o,t) + ('Uu‘p(0)) + ((g,‘p))(o,t) 

((ut,‘Pt))(o,t) - (ut(t),‘p(t)) - ((Vut,V(p))(o,t) - ((Vu, V(/5))(o,t) 

+ A ((«,(/?)) (o,t) + (wi,<p(0)) + ((g,<p))(o,t)- (3.16) 


A crucial point in our argument is that the left hand side tends, with no need of extracting a further 
subsequence, to a linear and continuous functional on Vt that acts on ip as specified by the right hand 
side. Noting as such a functional, we have in other words 

weakly in Vj. (3.17) 

Moreover, (|3.16l) can be restated as 

- ((Mt,<Pt))(o.t) + {u{t),ip{t)) + ((Vut, V(/5))(o,t) + ((Vu, V(p))(o,t) 

+ ({^(t),7’))(o,t) - A((M,(/?))(o,t) = (mi, 7’(0)) + ((5,7’))(o,t)- (3-18) 

Hence, (I2.31|l and (|2.30|1 . which is a particular case of it, are proved. Note now that, from (13.Mali , it 
also follows 

T[t) weakly star in X^. (3.19) 

and also this convergence holds with no need of extracting further subsequences. Indeed, the limit of 
the whole (sub)sequence is already identified as ^(t) on the dense subspace At fl V*. This also implies 
that the measure T[t) represents ^(t) on Xt in the sense of (12.1811 . Using the fact that for any ip G Vt,o 
the extension tp lies in V, it is easy to check that the functionals ^ and ^(t) are “compatible”. Hence, 
we have checked points (a) and (bl)-(b2) of Definition 12.21 of weak solution. 

Let us now show relation (I2.35I1 . i.e., point (e) of Definition 12.21 To this aim, we write (12.3111 
with if = u for s,t G (0,T] and take the difference. Note that the choice tp = u is admissible since 
u GV. We then infer 


- + (“t(^):“(^)) + / (VMt,Vu) dr + 

J S 

+ ((?(t)>w))(o,t) - ((^(s),w))(o,s) - = (ut(s),M(s)) + / / gu dx At. (3.20) 

J s J Q 


Then, computing explicitly the integral on the left hand side, (|2.35ll readily follows. 

Step 4. Identification of To conclude our proof we need to identify ^ (and ^(t)) in the sense of 
the weak constraint (12.3311 . This will give (c) of Definition l2.21 We start working on and, to get the 
identification, we shall implement the so-called Minty’s trick in the duality between V' and V. This 
corresponds to checking the following two conditions: 


(i) There holds the lim sup-inequality 

linisiip((/3^(u^),u^)} < U,u)). 

e\0 


(3.21) 


(ii) The operators suitably converge to in such a way that (|2.33ll may follow as a consequence 

of (ixim . 

We start by checking property (i), postponing the discussion regarding the correct notion of conver¬ 
gence for (ii) and its implications. Writing (1^ for ip = and t = T, we obtain 

{{P^{u%u^)) = ||<||i2(o.T;ff) - {ut{T),u^T)) + (nf,u§) - i||Vu^(r)|p + i||Vu^f 

~ + ((Si'*^^))- (3.22) 
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Now, thanks to (I3.15a|) . (I3.15cll . (I3.15dl) . p.l5g| , and (I3.15ell . we see that the limsup (as e \ 0 along 
a proper subsequence) of the right hand side is less or equal than 

- iMT),u{T)) + (mi,uo) - ^||Vu(r)f + i||Vuof - ||Vu||^ 2 ( 0 ,T;ff) 

+ M\u\\‘l2{o,t-,h) (3.23) 

Hence, using (12.351) written for t = T and s = 0, we see that the above expression is equal to {{^,u)). 
Therefore, (|3.2ip is proved. 

Let us now switch to discussing (ii), which requires the introduction of some additional ma¬ 
chinery. We present it by following the lines of the book by Attouch [1]. At first, we observe that the 
restriction to V of the function /3^ can be seen as a monotone operator from V to V' (once one works 
in the parabolic Hilbert triplet V C H C V')- Indeed, if u S V, then f5^{v) € H C V' by the Lipschitz 
continuity of /3^. Hence, for any ui,U 2 G V, we have 

iWM - (3 ^{vi),V 2 - vi)) = ((/ 3 ^(u 2 ) - l3‘^{vi),V 2 - vi)) > 0. (3.24) 

Moreover, if u, z G V, then, by definition of subdifferential, 

{{l3^v\z-v)) = {iP^v),z-v))< r [ {f{z)-f{v))dxdt = mz)-mv). (3.25) 

Jo Jn 

In other words, we have the graph inclusion 

/3^|v C (3.26) 

where the notation used on the right hand side stands for the subdifferential of J^\\! with respect to 
the duality pairing between V' and V. By the standard theory of subdifferentials, this is a maximal 
monotone operator from V to 2 ^ , which includes (in the sense of graphs) the (monotone, but not 
necessarily maximal) operator /3'^Iv¬ 
in view of the fact that the family of functionals (defined on V and taking values 

in [0,-l-oo)) is increasing as e decreases to 0, applying [H Thm. 3.20], we obtain 

J"|v^ sup J^\v (3.27) 

£ 6 ( 0 , 1 ) 

in the sense of Mosco convergence (that is Gamma-convergence both in the strong and in the weak 
topology of V). Moreover, by the monotone convergence theorem it is readily seen that the functional 
on the right hand side coincides in fact with J\\;- Hence, owing to [U Thm. 3.66], the family of 
maximal monotone operators identified with the family of their graphs in the product 

space V X V', converges in the sense of graphs (cf. [TJ Def. 3.58]) to dvy'J\v = Pw Namely, 

y[x;y] G dvy'J\v, 3[x^;?/'^] G dvy'J‘^\v such that [a;®;y^] [x]y\ strongly in V x V. (3.28) 

Hence, in view of the facts that [m®;/3®(u^)] G J^\v (thanks to (13.261) 1. [u^;,8®(m®)] ^ [u;^] weakly 
in V X V' (thanks to (|3.15all and p.l5j[ )), and to the limsup-inequality (13.211) . we may apply [TJ 
Prop. 3.59], yielding that [u;^] G dvyij^\v = l^w Hence, (12.331) is proved. 

To conclude this part, we need to prove that G To this aim, it is sufficient to 

adapt the above argument by working on the subinterval (0,t). Indeed, relations (13.91) and (12.351) (the 
latter for s = 0) hold on any subinterval (0,t). Moreover, we can take advantage of (I3.15al) (whose 
analogue obviously holds also on subintervals) and (13.171) . 

Step 5. Further properties of solutions. Let us start proving that inequality (13.211) is in fact an 
equality. Indeed, owing to (13.281) . there exist [x’^\y‘^] G dvy'J‘^\v such that [x’^\y‘^] [u;^] strongly 

in V X V'. Hence, noting that, by monotonicity, 

0<((/3^(u^)-^/^u^-a;^)), (3.29) 
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taking the liminf as e \ 0, and recalling (I3.21L we obtain 


lim ((/3®(u®), u®)) = U, u)). (3.30) 

e\0 

As a consequence, the limit of the right hand side of (13.221) exists and coincides with (13.2311 . In view 
of the fact that convergence of most terms of (13.2211 is already known from the previous estimates we 
get in particular that 


lim Q||Vw^(T)f + = i||Vw(r)||2 + (3.31) 

As before, this argument can be repeated on any subinterval (0,t). Hence, recalling (I3.15all and 
(13.1 Sell) , we finally arrive at 


—>■ u strongly in A^(0, T; H), (3.32) 

u^{t) — u{t) strongly in V for all t £ [0,T]. (3.33) 


Next, let us show that, under assumption (12.4211 . the additional regularity (12.4311 holds. To this 
aim, we go back to the approximate problem, and, in the spirit of [28], we test (13.611 by Au^. Indeed, 
u® has sufficient smoothness in order for this procedure to be admissible (cf. (13.711 1. Integrating by 
parts, and using the monotonicity of /3®, we then easily infer 

^{ul,Au‘^) + = (Am^ +g,Au^) + (uf,AMf). (3.34) 

By some further integration by parts and using Holder’s and Young’s inequalities (and the definition 
of the operator A), the right hand side can be easily estimated as follows: 


{Xu^+g,Au^) + {ulAul)<\\\Au^\^ + C{\\gr + \\u^\^) + \\Wuir. (3.35) 

Here and below, C > 0 is a constant independent of e. Hence, integrating (13.341) over (0,t) for 
arbitrary t G (0,T], using (I3.35|) . and recalling estimate (|3.12all . we easily obtain 


2{ul{t),Au'^{t)) + II 


+ / IIAm'II" ds^C+IIAu^l 


(3.36) 


Now, one can immediately check that, under assumption (12.4211 . if Uq is defined as in (m, then the 
right hand side of (13.361) is bounded independently of e. Hence, noticing that the left hand side is 
larger or equal than 

l\\Au^t)r - c\\umr + £ ( 3 . 37 ) 

where the second term is uniformly controlled due to (I3.12all . we readily arrive at 




(3.38) 


Letting e 0, we then infer 

u£ L°°{0,T-D{A)) (3.39) 

thanks to semicontinuity of norms with respect to weak convergence. Finally, (I2.43L i.e., weak con¬ 
tinuity of u with values in D{A), follows by combining (I3.39|) with the regularity u £ C([0,T];17) 
following from (12.2811 . and applying standard results. 

Eventually, we show that weak solutions constructed as limit points of {u®} also satisfy a form 
of the energy inequality. We start by proving it on intervals of the form [0,t], t £ (0,T]. To this aim, 
we write relation (|3.8I1 for s = 0 and take the liminf as e 0. Then, using standard semicontinuity 
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arguments together with relations (I3.2|) - (IH.3I1 . (I3.15all . and (Id.l5cl) - (l3.15d|) . it is not difficult to infer, 
for every t € (0, T], 

dx - ^||u(t)|p + i||Vu(t)|p + ||VMt||i 2 ( 0 .t;H) 

- + “ ^llwof + ^llVuof + ^ {g,Ut) dr. (3.40) 

Note in particular that relation 

J{u{t)) = [ j(M(t)) dec < liminf [ j^{u^{t)) dx = \imini J^{u^{t)) (3-41) 

Jn eNO Jq e\,0 

is a consequence of (I3.15cll and of the fact that the functionals J® converge to J in the sense of Mosco 
(cf. dJ Par. 3.3]) in the space H. Recalling (12.251) . (13.401) reduces to (12.411) in the case s = 0. Let us 
now consider a generic interval [s, t] for 0 < s < t < T and let us go back to (13.81) written for this 
choice of s,t. Let us take once more the liminf as £ \ 0. Then, the left hand side is treated exactly 
as before. On the other hand, when looking at the right hand side, it is easy to check that 

-■^l|w^(s)f+ ^||Vu^(s)|P + ^ (g,Mf) dr-s--^|| m(s)|P + i||VM(s)||2+ ^ (c,,ut) dr (3.42) 

thanks in particular to (l3.15aL and (13.331) . Next, thanks to (I3.15cl) . and up to extracting a further 
subsequence of e \ 0 , we have 

(3.43) 

for almost every choice of s € (0,T). Next, we need to control the component of the energy related 
with the constraint. Namely, we would like to prove that, at least for a.e. s € (0,T), 


limsup / j^{u^{s)) dx< j{u{s)) dx 

e\,o Jn Jn 


(3.44) 


(hence, coupling the above with (13.411) written for t = s, we would get convergence of that term). We 
can start noticing that 


[ f{u^s)) dx= [ {f{u^s))-f{u{s))) dx+ [ f{u{s)) dx. 


Moreover, 


fiu{s))dx^ / j{u[s)) dx 


by the monotone convergence theorem. Now, by definition of subdifferential, we may write 

f {f (u^is)) - finis))) dx< [ P%u%s)){u’^is)-u{s)) dxn%s), 

Jn Jn 

and we have to discuss the behavior of the functions /r®. First, we observe that 

[ Ai^(s) ds = ((/?®(m^),m^-I t))-J. 0, 

Jo 

the latter property following from (I3.30|) and p.lSjD . Moreover, thanks to (I3.47L we have 
fis) > Jfu^is)) - Jffs)) > Jfu^s)) - J(m(s)), 


(3.45) 

(3.46) 

(3.47) 

(3.48) 

(3.49) 


whence liminf£\^o fif > 0 thanks to (13.411) . Hence, we have in particular limE\^o(Ai‘^) (s) = 0, (•)■ 
denoting the negative part. Moreover, from (I3.49|) . (13.411) and (l3.12cD we infer 


fis) > —J{u{s)) > — liminf (it®(s)) > —M, 
( 5'\0 


(3.50) 
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for all e G (0,1) and s € (0,T]. Hence, by the dominated convergence theorem we obtain that 
(/r®)“ —>■ 0 in L^(0,T). Consequently, thanks to (I3.48|) . we conclude that 0 in L^(0,T). Hence, 
up to a subsequence, /r® —>■ 0 almost everywhere in (0, T), whence (13.441) follows. This actually implies 
(|2.41l) for almost every s G (0,T) and every t G (s,T], as desired. The proof of Theorem 12.51 is 
concluded. 

Remark 3.2. If the source term g is 0, from (12.4111 follows in particular that the energy loss in the 
time interval (s,t) is at least as large as the dissipation term V{s,t) := ||Vut||| 2 (s Of course, 
as commonly occurs situations characterized by bad regularity, the energy dissipated may be in fact 
strictly larger than V(s,t). Indeed, we may observe that proving equality in (13.401) appears out of 
reach in the present regularity setting. 

Remark 3.3. In view of our strategy of proof for Theorem l2.51 we can give some further observation 
complementing Bemark 12.41 Hence, let us go back to the “toy problem” (12.3611 . for example with 
P = 3/[_i 1 ] (but our consideration also apply to different choices of /3). Then, implementing our 
regularization method we get the equation 

Wtt+/3''('«'') = 0. (3.51) 

Setting := uf, the shape of solution trajectories of the 2D ODE system associated to (13.511) in 
the phase space (m®,u®) can be easily described. In particular, since (in this spatially homogeneous 
setting) no dissipation occurs, trajectories are periodic. Moreover, we may notice that, for e 0, 
(u®, u®) converges in a suitable way to a couple (u, v), where v = Ut and u solves (j2.36ll . Clearly, (u, v) 
is also a periodic trajectory and its image in the phase space lies in some level set {j{u) + v‘^/2 = c}, 
c > 0, of the “energy” functional. In particular, whenever D{j) = [—1,1] (as happens in the case of the 
indicator function j = /[_i ij, and also for the “logarithmic potential” mentioned in the introduction), 
such level sets are (at least for large initial energy, i.e., for large values of c) not connected. Namely, 
their shape determines the jumps of Ut (which, consequently, cannot occur in an “arbitrary” way). 
Note also that taking different choices for the approximations of P does not modify the shape 
of {u,v). Of course it is clear that, in the case of our equation (HU, the situation is much more 
complicated than for (12.361) in view of the infinite-dimensional setting. However, the fact that our 
weak solutions u are still built as limit points of families u® solving a very natural regularization of 
the equation suggests that the jumps of Ut occurring in the limit may be in some sense “physical”, 
i.e., they are determined by the fact that uf, as e \ 0, may tend to develop discontinuities. In other 
words, the occurrence of “spurious” jumps of Ut (as are the somehow “arbitrary” jumps described 
in Remark 12. 4|) should be excluded in view of the fact that our weak solutions descend from the 
approximation scheme. 

Remark 3.4. Let us give some further observation complementing Remark 1 2.3 1 Again, we consider, 
just for simplicity, the “toy” model (12.361) : however, our considerations also apply to the original 
equation (HU. Actually, from our approximation argument we know that, for any t G (0,T], (a 
subsequence of) P^{u^) (weakly star) converges to a measure 7(t) on Qt (in particular, we have 
convergence to some T on the whole interval). In the toy case, of course, Qt = [0,t]; moreover, 
we are allowed to identify 7(t) = C(t) £ because Sobolev functions are continuous in ID. Let us 
now consider the particular case when P^{u^) is supported in some interval of the form \t^ — e^t^ -\- e] 
and is 0 outside that interval. Then, assuming that t^ converges to some point t G (0,r) as £ \ 0, 
and P^{u^) “spikes” around in a proper way, it may happen that P^{u^) (weakly star) converges to 
T = St (the Dirac delta concentrated in t) in X'. This kind of behavior may be (possibly) driven for 
instance by inserting a nonzero forcing term g in the equation. Then, in the case when, for instance, 
= t — 2e, it turns out that the singularity of T develops before t. Consequently, P^{u^) also converges 
to St in Xp In particular, (12.311) holds in [0,t] with Qt) = St- On the other hand, ii = t + 2e, i.e., 
the singularity of T develops after t, in that case P^{u^) converges to 7(t) = 0 in XQ whence (12.311) 
holds in [0,t] with Qt) = 0- Note that this happens in spite of the fact that the limit measure T over 
the whole [0,T] is the same in the two cases. This fact suggests that the formulation (12.311) on the 
subinterval (0, t) contains some additional information that cannot be simply inferred by restricting 
the global formulation (12.301) . This is the reason why we decided to include (b2) in our existence 
theorem. 
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Remark 3.5. Let t be one of the (at most countably many) jump points of Ut- Then, both the point 
value of Ut at t and the occurrence of concentration phenomena for the measure 7(t) at the same point 
also depend on the choice of the approximating problem (i.e., of actually our argument works 
provided that (3^ is smooth and converges to /3 in the graph sense) and of the selection of converging 
subsequences via Helly’s theorem. This can be seen again by looking at the “toy equation” (12.3611 
with P = and initial values uq = 0 and ui = 1. Then, we know that the (first) jump of Ut 

occurs at t = 1. Let us now consider the approximation (13.5111 with the choice 


{ 0 if |r| < r^, 

e-2(r-r^) if r > 

e-2(r + r^) if r < 


(3.52) 


where, for any e S (0,1), one may choose (in an arbitrary way) in the interval [1 — evr, 1]. It is then 
clear that, whatever are the chosen values of r^, (3^ tends to P = 9/[_i ij in the sense of graphs as 
e \j 0. Hence, our limit problem is the desired one. Let us notice that, for r > 0, we have 


In particular, j^(l) = 0 if = 1, whereas 

V n o 

(o, 


(1 - r '^)2 

/(I) = e 


2 e 2 


TT 

T 


if G [1 — ETT, 1). 


(3.53) 


(3.54) 


Then, it®(l) = 1 if = 1, whereas in case r® G [1 — ett, 1) one can easily compute u^{t) = r® + 
Esin for t G (r®,r® + ett], whence ul{t) = cos and uf{l) = cos Hence, choosing 

appropriately (and somehow “wildly”) in the interval [1 — ett, 1] as e varies in (0,1), one may obtain 
the effect that for any number i G [—1,1] there exists a subsequence e^ \ 0 such that M("(1) tends 
to i. The use of Helly’s theorem selects one of these subsequences and determines the limit value 
Ut(l) = i (and, in turn, how the limit measures 7(s) concentrate at the jump point t = 1). 


Remark 3.6. As observed in the previous Remark, is not represented, in general, by the restricion 
of the measure T to the set Qt- However, we can give a more explicit characterization of this restriction 
in the following sense. From (I3.15h|) we have Ut G BV{0,T; X), where we may assume X be the dual 
of a separable space. Hence, for all times t G [0,T] there exists (in the weak star topology of X) the 
limit 


Ut{t^) := rc*—lim Ut{s). 


t+ 


Moreover this value coincides with the weak star limit 

1 /■'* 

w*—lim 


s —S — t 


Ut{r) dr. 


(3.55) 


(3.56) 


In particular the limits above must hold with respect to the weak topology of H, since Ut is bounded 
in H uniformly in time. Let us now write (12.301) with ip replaced by phs G A, with G A (~i V and 
hs ; [0,r] — )■ [0, 1 ] be the function such that hg = 1 on [0,t], hg = 0 on [s,T], and hg be affine in [t,s]. 
We obtain 


- ((wt,¥>t))(o.t) - iut,pthg))(t,s) 


s — t 


iut,p))(t,s) + ((Vut, V((/j/ig))) + {{Vu,V{phg))) 


dT - X{{u,phg)) = {ui,p{0)) + {{g,phg)). 


Letting s \t, we see that the third term tends to while the other terms pass to the 

limit thanks to the dominated convergence theorem and the fact that phg —>■ pX]fj t] pointwise, so in 
particular T-almost everywhere. We then obtain 


- {(.ut,<Pt)){o,t) + (ut(t+), (/?(<)) + ((Vut, V<p))(0,t) + ((Vu, V(/5))(0,t) 


+ 



p d(rX[o,t]) - A((u,v3))(o,t) 


(m 1,(/3(0)) + ((g,(/7))(o,t). 
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Comparing with ( 12 .Hill , we deduce that ^(t) is represented by the restriction of T to the closed set Qt 
whenever the pointwise value Ut{t) coincides with which happens in fact in the complementary 

of a countable set of times. In other words, in that case we have 'Tl.-^= Tit)- 

Remark 3.7. Relation (12.4111 implies in particular that, at least when g = 0, the energy functional 
coincides almost everywhere with a nonincreasing function. In a sense this fact provides an additional 
criterion for selecting which are the “admissible” jumps of Ut (cf. Remark 13.31) . Namely, jumps may 
occur only in such a way that they do not increase the total energy of the system. For 5 7 ^ 0 similar 
considerations hold, up to the fact that g acts somehow as an additional energy source. 

Remark 3.8. It is maybe also worth stressing that Theorem 12.51 states the existence of at least one 
weak solution satisfying the properties detailed above. Due to nonuniqueness, there may well exist 
“spurious” solutions having worse properties. For example they may be constructed in such a way 
that the time derivative Ut admits somehow “nonphysical” jumps. However our procedure shows 
that every weak solution that is a limit point of our natural regularization scheme is “physical” (for 
example, in view of (I2.41L energy-increasing jumps cannot occur). 
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